Analysis of cancer-associated fibroblasts related genes identifies COL11A1 associated with lung adenocarcinoma prognosis.

BACKGROUND: The treatment of lung adenocarcinoma is difficult due to the limited therapeutic options. Cancer-associated fibroblasts play an important role in the development of cancers. This study aimed to identify a promising molecular target associated with cancer-associated fibroblasts for the treatment of lung adenocarcinoma. METHODS: The Cancer Genome Atlas lung adenocarcinoma dataset was used to screen hub genes associated with cancer-associated fibroblasts via the EPIC algorithm and Weighted Gene Co-expression Network Analysis. Multiple databases were used together with our data to verify the differential expression and survival of COL11A1. Functional enrichment analysis and the single-cell TISCH database were used to elucidate the mechanisms underlying COL11A1 expression. The correlation between COL11A1 and immune checkpoint genes in human cancers was also evaluated. RESULTS: Using the EPIC algorithm and Weighted Gene Co-expression Network Analysis, 13 hub genes associated with cancer-associated fibroblasts in lung adenocarcinoma were screened. Using the GEPIA database, Kaplan-Meier Plotter database, GSE72094, GSE75037, GSE32863, and our immunohistochemistry experiment data, we confirmed that COL11A1 overexpresses in lung adenocarcinoma and that high expression of COL11A1 is associated with a poor prognosis. COL11A1 has a genetic alteration frequency of 22% in patients with lung adenocarcinoma. COL11A1 is involved in the extracellular matrix activities of lung adenocarcinoma. Using the TISCH database, we found that COL11A1 is mainly expressed by cancer-associated fibroblasts in the tumor microenvironment rather than by lung adenocarcinoma cells. Finally, we found that COL11A1 is positively correlated with HAVCR2(TIM3), CD274 (PD-L1), CTLA4, and LAG3 in lung adenocarcinoma. CONCLUSION: COL11A1 may be expressed and secreted by cancer-associated fibroblasts, and a high expression of COL11A1 may result in T cell exhaustion in the tumor microenvironment of lung adenocarcinoma. COL11A1 may serve as an attractive biomarker to provide new insights into cancer therapeutics.

Correction: Zheng et al. Glycolysis-Related SLC2A1 Is a Potential Pan-Cancer Biomarker for Prognosis and Immunotherapy. Cancers 2022, 14, 5344.

The authors wish to make the following corrections to this paper [...].

Addition of bevacizumab to EGFR tyrosine kinase inhibitors in advanced NSCLC: an updated systematic review and meta-analysis.

Background: The synergistic effects of antiangiogenic inhibitor bevacizumab and epidermal growth factor receptor-tyrosine kinase inhibitors (EGFR-TKI) therapy were encouraging in patients with EGFR-mutant advanced NSCLC, though some controversy remains. The specific subgroup of patients who might benefit most from the EGFR-TKI and bevacizumab combination therapy is yet to be determined. Methods: Randomized clinical trials (RCTs) that had compared the clinical efficacy of EGFR-TKI and bevacizumab combination therapy with EGFR-TKI monotherapy in treating EGFR-mutant advanced NSCLC patients published before 23 December 2022 were searched in the Cochrane, PubMed and Embase. We performed a meta-analysis for the overall survival (OS), progression-free survival (PFS), objective response rate (ORR), and treatment-related adverse events with a grade equal or more than 3 (grade≥3 TRAEs). Subgroup analyses of PFS and OS stratified by clinical characteristics and treatment were conducted. Results: We included 10 RCTs involving 1520 patients. Compared with EGFR-TKI monotherapy, addition of bevacizumab to EGFR-TKI resulted in a significantly higher PFS (hazard ratio (HR) = 0.74, 95% confidence interval (95% CI): 0.62-0.87)) and ORR (risk ratio (RR) = 1.07, 95% CI: 1.01-1.13). However, no significant difference in OS (HR = 0.96, 95% CI: 0.83-1.12) was noticed. Patients with EGFR-mutant advanced NSCLC receiving combination therapy showed PFS improvement regardless of gender (male or female), Eastern Cooperative Oncology Group performance status (0 or 1), baseline central nervous system (CNS) metastasis (presence or absence) and EGFR mutation type (19del or 21L858R). Subgroup analyses showed that, with the treatment of bevacizumab and EGFR-TKI, patients who ever smoked achieved significantly better OS and PFS benefits (HR = 0.68, 95% CI: 0.48-0.95; HR = 0.59, 95% CI: 0.46-0.74, respectively), and those aged <75 years and the Asian population had significantly prolonged PFS (HR = 0.69, 95% CI: 0.52-0.91; HR = 0.71, 95% CI: 0.58-0.87; respectively). The superiority of EGFR-TKI and bevacizumab combination therapy against EGFR-TKI monotherapy in improving PFS was more significant in the erlotinib regimen subgroup. The risk of grade≥3 TRAEs was remarkably higher in the combination therapy group (HR = 1.73, 95% CI: 1.39-2.16). Conclusion: Addition of bevacizumab to EGFR-TKI therapy provided significantly better PFS and ORR for EGFR-mutant advanced NSCLC patients, though with higher risk of grade≥3 TRAEs. Patients who ever smoked, aged <75 years, and Asian population might benefit more from the combination regimen. Systematic Review Registration: This systematic review and meta-analysis was registered in the PROSPERO database (CRD42023401926).

Glycolysis-Related SLC2A1 Is a Potential Pan-Cancer Biomarker for Prognosis and Immunotherapy.

SLC2A1 plays a pivotal role in cancer glycometabolism. SLC2A1 has been proposed as a putative driver gene in various cancers. However, a pan-cancer analysis of SLC2A1 has not yet been performed. In this study, we explored the expression and prognosis of SLC2A1 in pan-cancer across multiple databases. We conducted genetic alteration, epigenetic, and functional enrichment analyses of SLC2A. We calculated the correlation between SLC2A1 and tumor microenvironment using the TCGA pan-cancer dataset. We observed high expression levels of SLC2A1 with poor prognosis in most cancers. The overall genetic alteration frequency of SLC2A1 was 1.8% in pan-cancer, and the SLC2A1 promoter was hypomethylation in several cancers. Most m6A-methylation-related genes positively correlated with the expression of SLC2A1 in 33 TCGA cancers. Moreover, SLC2A1 was mainly related to the functions including epithelial-mesenchymal transition, glycolysis, hypoxia, cell-cycle regulation, and DNA repair. Finally, SLC2A1 positively associated with neutrophils and cancer-associated fibroblasts in the tumor microenvironment of most cancers and significantly correlated with TMB and MSI in various cancers. Notably, SLC2A1 was remarkably positively correlated with PD-L1 and CTLA4 in most cancers. SLC2A1 might serve as an attractive pan-cancer biomarker for providing new insights into cancer therapeutics.

Mining the Prognostic Role of DNA Methylation Heterogeneity in Lung Adenocarcinoma.

Purpose: DNA methylation heterogeneity is a type of tumor heterogeneity in the tumor microenvironment, but studies on the identification of the molecular heterogeneity of the lung adenocarcinoma genome with respect to DNA methylation sites and their roles in lung cancer progression and prognosis are scarce. Methods: Prognosis-associated DNA methylation subtypes were filtered by the Cox proportional hazards model and then established by unsupervised cluster analysis. Association analysis of these subtypes with clinical features and functional analysis of annotated genes potentially affected by methylation sites were performed. The robustness of the model was further tested by a Bayesian network classifier. Results: Over 7 thousand methylation sites were associated with lung adenocarcinoma prognosis. We identified seven molecular methylation subtypes, including 630 methylation sites. The subtypes yielded the most stable results for differentiating methylation profiles, prognosis, and gene expression patterns. The annotated genes potentially affected by these methylation sites are enriched in biological processes such as morphogenesis and cell adhesion, but their individual impact on the tumor microenvironment and prognosis is multifaceted. Discussion. We revealed that DNA methylation heterogeneity could be clustered and associated with the clinical features and prognosis of lung adenocarcinoma, which could lead to the development of a novel molecular tool for clinical evaluation.

Establishing a TNM-like risk classification for metachronous second pulmonary adenocarcinoma in patients with previously resected pulmonary adenocarcinoma.

Background: For metachronous second pulmonary adenocarcinoma (msPAD) in patients with resected PAD, the method to distinguish tumour clonality has not yet been well established, which makes it difficult to determine accurate staging and predict prognosis. Methods: Patients received surgery for the primary and encountered msPAD were recruited into the Surveillance, Epidemiology, and End Results database. We extracted overall survival 1 (OS1) for the primary, overall survival 2 (OS2) for the msPAD, and defined interval survival as the interval time between the first and second PAD. Based on the nomogram and recursive partitioning analysis, a tumor, node, metastasis staging system (TNM)-like risk stratification system was established for OS2 on the premise of suspending the dispute of tumor clonality. Results: A total of 1,045 patients were identified. There is no significant association between interval survival and OS2. A TNM-like risk stratification system was established based on the independent pathological factors for prognosis, including tumor diameter (2nd), node metastasis (2nd), grade (2nd), and extrapulmonary metastasis (2nd). The proposed risk stratification system present well capacity in predicting and stratifying prognosis. Compared with the TNM stage system, the proposed risk stratification system presents a smaller Akaike information criterion (AIC) but larger c-index, and generates higher accuracy to predict prognosis at 160 months of follow-up according to the time-dependent receiver operating curve (ROC) curve. Conclusions: In conclusion, the TNM-like risk stratification appears to be suitable for prognostic prediction and risk stratification for msPAD patients with former PAD resection. This model validates and refines the known classification rules based on the easily collected variables, and highlights potentially clinical implications.

Establishing a prognostic model for metachronous second squamous cell lung cancer in patients with resected squamous cell lung cancer.

BACKGROUND: For metachronous second pulmonary squamous cell carcinoma (msPSC) in patients with resected PSC, the method to distinguish tumour clonality has not yet been well established, which makes it difficult to determine accurate staging and predict prognosis. METHODS: Patients who underwent surgery for first PSC and encountered msPSC were recruited from the Surveillance, Epidemiology, and End Results (SEER) database. We extracted overall survival 1 (OS1) for the first PSC, overall survival 2 (OS2) for msPSC, and interval survival for the time interval between the first and second PSC. The nomogram was calibrated for OS2, and recursive partitioning analysis (RPA) was performed for risk stratification. RESULTS: A total of 617 patients were identified. Several independent prognostic factors were identified and integrated into the nomogram for OS2, including gender, age (2nd), nodal status (1st), node metastasis (2nd), and extrapulmonary metastasis (2nd). The calibration curves showed optimal agreement between the predictions and actual observations, and the c-index was 0.678. Surgery was associated with longer survival for msPSC patients. The prognosis of sublobectomy was comparable and inferior to that of lobectomy in the low- and moderate-risk groups, respectively. Radiotherapy was associated with better outcomes in patients who did not undergo surgery. CONCLUSIONS: The RPA-based clinical nomogram appears to be suitable for the prognostic prediction and risk stratification of OS2 in msPSC. This practical system may help clinicians make decisions and design clinical studies.

Systematic analysis of myocardial immune progression in septic cardiomyopathy: Immune-related mechanisms in septic cardiomyopathy.

Background: The immune infiltration and molecular mechanisms underlying septic cardiomyopathy (SC) have not been completely elucidated. This study aimed to identify key genes related to SC and elucidate the potential molecular mechanisms. Methods: The weighted correlation network analysis (WGCNA), linear models for microarray analysis (LIMMA), protein-protein interaction (PPI) network, CIBERSORT, Kyoto Encyclopedia of Genes and Genomes pathway (KEGG), and gene set enrichment analysis (GSEA) were applied to assess the key pathway and hub genes involved in SC. Results: We identified 10 hub genes, namely, LRG1, LCN2, PTX3, E LANE, TCN1, CLEC4D, FPR2, MCEMP1, CEACAM8, and CD177. Furthermore, we used GSEA for all genes and online tools to explore the function of the hub genes. Finally, we took the intersection between differential expression genes (DEGs) and hub genes to identify LCN2 and PTX3 as key genes. We found that immune-related pathways played vital roles in SC. LCN2 and PTX3 were key genes in SC progression, which mainly showed an anti-inflammatory effect. The significant immune cells in cardiomyocytes of SC were neutrophils and M2 macrophages. Conclusion: These cells may have the potential to be prognostic and therapeutic targets in the clinical management of SC. Excessive anti-inflammatory function and neutrophil infiltration are probably the primary causes of SC.

Multi-factor regulatory network and different clusters in hypertrophic obstructive cardiomyopathy.

BACKGROUND: Practical biosignatures and thorough understanding of regulatory processes of hypertrophic obstructive cardiomyopathy (HOCM) are still lacking. METHODS: Firstly, public data from GSE36961 and GSE89714 datasets of Gene Expression Omnibus (GEO), Gene database of NCBI (National Center of Biotechnology Information) and Online Mendelian Inheritance in Man (OMIM) database were merged into a candidate gene set of HOCM. Secondly, weighted gene co-expression network analysis (WGCNA) for the candidate gene set was carried out to determine premier co-expressed genes. Thirdly, significant regulators were found out by virtue of a multi-factor regulatory network of long non-coding RNAs (lncRNAs), messenger RNAs (mRNAs), microRNAs (miRNAs) and transcription factors (TFs) with molecule interreactions from starBase v2.0 database and TRRUST v2 database. Ultimately, HOCM unsupervised clustering and "tsne" dimensionality reduction was employed to gain hub genes, whose classification performance was evaluated by a multinomial model of lasso logistic regression analysis binded with receiver operating characteristic (ROC) curve. RESULTS: Two HOCM remarkably-interrelated modules were from WGCNA, followed by the recognition of 32 crucial co-expressed genes. The multi-factor regulatory network disclosed 7 primary regulatory agents, containing lncRNAs (XIST, MALAT1, and H19), TFs (SPI1 and SP1) and miRNAs (hsa-miR-29b-39 and has-miR-29a-3p). Four clusters of HOCM and 4 hub genes (COMP, FMOD, AEBP1 and SULF1) significantly expressing in preceding four subtypes were obtained, while ROC curve demonstrated satisfactory performance of clustering and 4 genes. CONCLUSIONS: Our consequences furnish valuable resource which may bring about prospective mechanistic and therapeutic anatomization in HOCM.

Asporin is a Potential Promising Biomarker for Common Heart Failure.

Heart failure (HF) is the end-stage of various diseases, especially ischemic cardiomyopathy (ICM) and dilated cardiomyopathy (DCM). We aimed to investigate the common molecular mechanism of ICM and DCM. Differentially expressed genes (DEGs) of ICM or DCM samples compared with control were identified in GSE1869, GSE5406, GSE57338, GSE79962, GSE116250, and GSE46224 datasets. Functional enrichment analysis and protein-protein network analysis of the coregulated DEGs in at least four datasets were performed using the online tools of DAVID, the Metascape database, and the STRING database. Hub genes of HF were identified and validated by western blotting (WB) and immunohistochemistry in our tissue microarray (TMA). Seventy-four coregulated ICM and 126 coregulated DCM relevant DEGs were identified. Moreover, 59 common genes between ICM and DCM relevant DEGs were obtained, which were mainly involved in cardiac fibrosis and several signal pathways, such as Wnt signal pathway, PI3K-Akt signal pathway, and HIF-1A signal pathway. Among the six hub genes with top degrees, asporin (ASPN) had a relatively higher correlation with LVEF. Finally, TMA and WB results revealed that the ASPN protein was significantly increased in ICM and DCM left ventricular samples. The present study revealed some common molecular mechanisms of HF with different causes. Furthermore, ASPN may be a potential promising biomarker for HF.

FLNC and MYLK2 Gene Mutations in a Chinese Family with Different Phenotypes of Cardiomyopathy.

Mutations in the sarcomeric protein filamin C (FLNC) gene have been linked to hypertrophic cardiomyopathy (HCM), as they have been determined to increase the risk of ventricular arrhythmia and sudden death. Thus, in this study, we identified a novel missense mutation of FLNC in a Chinese family with HCM, and, interestingly, a second novel truncating mutation of MYLK2 was discobered in one family member with different phenotype.We performed whole-exome sequencing in a Chinese family with HCM of unknown cause. To determine and confirm the function of a novel mutation of FLNC, we introduced the mutant and wild-type gene into AC16 cells (human cardiomyocytes): we then used western blotting to analyze the expression of FLNC in subcellular fractions, and confocal microscope to observe the subcellular distribution of the protein.As per our findings, we were able to identify a novel missense single nucleotide variant (FLNC c.G5935A [p.A1979T]) in the family, which segregates with the disease. FLNC expression levels were observed to be equivalent in both wild-type and p.A1979T cardiomyocytes. However, the expression of the mutant protein has resulted in cytoplasmic protein aggregations, in contrast to wild-type FLNC, which was distributed in the cytoplasm and did not form aggregates. Unexpectedly, a second truncating mutation, NM_033118:exon8:c.G1138T:p.E380X of the MYLK2 gene, was identified in the mother of the proband with dilated cardiomyopathy, which was not found in other subjects.We then identified the FLNC A1979T mutation as a novel pathogenic variant associated with HCM in a Chinese family as well as a second causal mutation in a family member with a distinct phenotype. The possibility that there is more than one causal mutation in cardiomyopathy warrants clinical attention, especially for patients with atypical clinical features.

Systematic analysis of molecular mechanisms of heart failure through the pathway and network-based approach.

AIMS: The molecular networks and pathways involved in heart failure (HF) are still largely unknown. The present study aimed to systematically investigate the genes associated with HF, comprehensively explore their interactions and functions, and identify possible regulatory networks involved in HF. MAIN METHODS: The weighted gene coexpression network analysis (WGCNA), crosstalk analysis, and Pivot analysis were used to identify gene connections, interaction networks, and molecular regulatory mechanisms. Functional analysis and protein-protein interaction (PPI) were performed using DAVID and STRING databases. Gene set variation analysis (GSVA) and receiver operating characteristic (ROC) curve analysis were also performed to evaluate the relationship of the hub genes with HF. KEY FINDINGS: A total of 5968 HF-related genes were obtained to construct the co-expression networks, and 18 relatively independent and closely linked modules were identified. Pivot analysis suggested that four transcription factors and five noncoding RNAs were involved in regulating the process of HF. The genes in the module with the highest positive correlation to HF was mainly enriched in cardiac remodeling and response to stress. Five upregulated hub genes (ASPN, FMOD, NT5E, LUM, and OGN) were identified and validated. Furthermore, the GSVA scores of the five hub genes for HF had a relatively high areas under the curve (AUC). SIGNIFICANCE: The results of this study revealed specific molecular networks and their potential regulatory mechanisms involved in HF. These may provide new insight into understanding the mechanisms underlying HF and help to identify more effective therapeutic targets for HF.

Global gene expression analysis using RNA-seq reveals the new roles of Panax notoginseng Saponins in ischemic cardiomyocytes.

ETHNOPHARMACOLOGICAL RELEVANCE: Panax notoginseng saponins (PNS), the main active ingredients of Panax notoginseng (Burkill) F.H.Chen, have been clinically used for cardiovascular diseases treatment in China as the Traditional Chinese Medicine (TCM) (Duan et al., 2017). Evidence demonstrated that PNS protected cardiomyocytes from myocardial ischemia, but the more underlying molecular mechanisms of the protective effect are still unclear. The aims of this study are to systematically know the function of PNS and discover new roles of PNS in ischemic cardiomyocytes. MATERIALS AND METHODS: To confirm PNS function on ischemic cardiomyopathy, we established in vitro myocardial ischemia model on H9C2 cardiomyocyte line, which was induced by oxygen-glucose depletion (OGD). Then RNA-seq was carried out to systematically analyze global gene expression. This study was aimed to systematically investigate the protective effect and more potential molecular mechanisms of PNS on H9C2 cardiomyocytes in vitro through whole-transcriptome analysis with total RNA sequencing (RNA-Seq). RESULTS: PNS exhibited anti-apoptotic effect in H9C2 cardiomyocytes in OGD-induced myocardial ischemia model. Through RNA-seq, we found that OGD affected expression profiling of many genes, including upregulated and downregulated genes. PNS inhibited cardiomyocyte apoptosis and death through rescuing cell cycle arrest, the DNA double-strand breakage repair process and chromosome segregation. Interestingly, for the canonical signaling pathways regulation, RNA-seq showed PNS could inhibit cardiac hypertrophy, MAPK signaling pathway, and re-activate PI3K/AKT and AMPK signaling pathways. Experimental data also confirmed the PNS could protect cardiomyocytes from OGD-induced apoptosis through activating PI3K/AKT and AMPK signaling pathways. Moreover, RNA-seq demonstrated that the expression levels of many non-coding RNAs, such as miRNAs and lncRNAs, were significantly affected after PNS treatment, suggesting that PNS could protect cardiomyocytes through regulating non-coding RNAs. CONCLUSION: RNA-seq systematically revealed different novel roles of Panax Notoginseng Saponins (PNS) in protecting cardiomyocytes from apoptosis, induced by myocardial ischemia, through rescuing cell cycle arrest and cardiac hypertrophy, re-activating the DNA double-strand breakage repair process, chromosome segregation, PI3K/Akt and AMPK signaling pathways and regulating non-coding RNAs.

Analysis of genes and underlying mechanisms involved in foam cells formation and atherosclerosis development.

BACKGROUND: Foam cells (FCs) play crucial roles in the process of all stages of atherosclerosis. Smooth muscle cells (SMCs) and macrophages are the major sources of FCs. This study aimed to identify the common molecular mechanism in these two types of FCs. METHODS: GSE28829, GSE43292, GSE68021, and GSE54666 were included to identify the differentially expressed genes (DEGs) associated with FCs derived from SMCs and macrophages. Gene Ontology biological process (GO-BP) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed by using the DAVID database. The co-regulated genes associated with the two origins of FCs were validated (GSE9874), and their expression in vulnerable atherosclerosis plaques (GSE120521 and GSE41571) was assessed. RESULTS: A total of 432 genes associated with FCs derived from SMCs (SMC-FCs) and 81 genes associated with FCs derived from macrophages (M-FCs) were identified, and they were mainly involved in lipid metabolism, inflammation, cell cycle/apoptosis. Furthermore, three co-regulated genes associated with FCs were identified: GLRX, RNF13, and ABCA1. These three common genes showed an increased tendency in unstable or ruptured plaques, although in some cases, no statistically significant difference was found. CONCLUSIONS: DEGs related to FCs derived from SMCs and macrophages have contributed to the understanding of the molecular mechanism underlying the formation of FCs and atherosclerosis. GLRX, RNF13, and ABCA1 might be potential targets for atherosclerosis treatment.

Role of the ADCY9 gene in cardiac abnormalities of the Rubinstein-Taybi syndrome.

BACKGROUND: Rubinstein-Taybi syndrome (RTS) is a rare, congenital, plurimalformative, and neurodevelopmental disorder. Previous studies have reported that large deletions contribute to more severe RTS phenotypes than those caused by CREBBP point mutations, suggesting a concurrent pathogenetic role of flanking genes, typical of contiguous gene syndromes, but the detailed genetics are unclear. RESULTS: This study presented a rare case of Rubinstein-Taybi (RT) syndrome with serious cardiac abnormalities. Based on the clinical and genetic analysis of the patient, the ADCY9 gene deletion was highlighted as a plausible explanation of cardiac abnormalities. In adcy9 morphant zebrafish, cardiac malformation was observed. Immunofluorescence study disclosed increased macrophage migration and cardiac apoptosis. RNA sequencing in zebrafish model highlighted the changes of a number of genes, including increased expression of the mmp9 gene which encodes a matrix metalloproteinase with the main function to degrade and remodel extracellular matrix. CONCLUSIONS: In this study, we identified a plausible new candidate gene ADCY9 of CHD through the clinical and genetic analysis of a rare case of Rubinstein-Taybi (RT) syndrome with serious cardiac abnormalities. By functional study of zebrafish, we demonstrated that deletion of adcy9 is the causation for the cardiac abnormalities. Cardiac apoptosis and increased expression of the MMP9 gene are involved in the pathogenesis.